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ABSTRACT 

Observations of narrow radio-emitting filaments near the Galactic center have been interpreted in previous 
studies as evidence of a pervasive vertical (i.e. _L to the Galactic plane) milliGauss magnetic field in the central 
^ 150 pc of the Galaxy. A simple cylindrically symmetric model for the equilibrium in this central region is 
proposed in which horizontal (i.e. || to the Galactic plane) magnetic fields embedded in an annular band of 
partially ionized molecular material of radius ~ 150 pc are wrapped around vertical magnetic fields threading low- 
density hot plasma. The central vertical magnetic field, which has a pressure that significantly exceeds the thermal 
pressure of the medium, is confined by the weight of the molecular material. The stability of this equilibrium 
is studied indirectly by analyzing a uniformly rotating cylinder of infinite extent along the z axis in cylindrical 
coordinates (r, 9,z), with low-density plasma and an axial magnetic field at r < 150 pc, high-density plasma 
and an azimuthal field at r > 150 pc, and a gravitational acceleration g* °^ r directed in the — f direction. Simple 
profiles are assumed for the density p, pressure p, and field strength B, with the sound speed and Alfven speeds <=<: r 
within the dense plasma. The density profile and gravity tend to destabilize the plasma, but the plasma tends to be 
stabilized by rotation and magnetic tension — since the interface between the high and low-density plasmas can not 
be perturbed without bending either the horizontal or vertical field. Normal modes proportional to e("«9+''^;z-'°') 
with k^ — and m^Q are studied. Such modes neither bend nor compress the axial field at r < 150 pc but allow 
compressions of the dense plasma along the azimuthal magnetic field that enhance the destabilizing role of gravity. 
It is shown analytically that when P = Snp/B^ is small and the dense plasma is supported against gravity primarily 
by rotation, the necessary and sufficient condition for stability to k^ = modes is \g\ < 2|i2|fl, where g = g* ~ ii^r 
is the effective gravity, Q. is the uniform angular velocity, and a is the sound speed in the dense plasma. Since the 
effective gravity is determined by the degree to which magnetic (and to a lesser degree pressure) forces support 
the dense plasma, the stability criterion gives an upper limit on the strength of the axial magnetic field, which is 
^ 1 mG for Galactic-center parameters. 



1. INTRODUCTION 

The central ~ 150 pc of the Galaxy host a phenomenon that 
is unique within the Galaxy: narrow, "vertical" (i.e. _L to the 
Galactic plane) filaments of radio emission that are typically a 
few tens of parsecs long and only a fraction of a parsec wide 
(Morris 1996). Polarization measurements indicate that the 
magnetic fields in the Radio Arc and Northern Thread are par- 
allel to these filaments (Tsuboi et al. 1986, Reich 1994, Tsuboi 
et al. 1995, Lang et al. 1999b), and this alignment is presumed 
to be typical of the other filaments as well. In virtually all cases 
the filaments appear to be in contact with a molecular cloud at 
some point along their length. The absence of bending of these 
filaments by the randomly moving clouds has been interpreted 
as evidence of a milliGauss lower limit to the field strength 
in the filaments (Yusef-Zadeh & Morris 1987a — this point is 
discussed further in section 7). The apparent impossibility of 
forming or confining narrow filaments with such tremendous 
magnetic pressures — well in excess of the thermal pressure of 
the medium — suggests that if a mG vertical field is present, it 
pervades the central ^ 150 pc of the Galaxy, and the filaments 
are simply those flux tubes that possess large populations of 
relativistic electrons (Morris & Serabyn 1996). 

A strong pervasive magnetic field at the Galactic center 
would be a natural consequence of radial inflow in the Galac- 
tic disk (Sofue & Fujimoto 1987, Morris & Serabyn 1996), and 
magnetic-pressure forces would over time cause the volume- 



filling magnetic flux tubes to rise buoyantly away from the 
Galactic mid-plane, forcing field lines that penetrate the mid- 
plane into a vertical orientation (Chandran et al. 2000). It is not 
clear, however, that such a scenario is inevitable, as the rate of 
radial inflow in the Galactic disk is not well known and since 
sufficient turbulent resistivity can inhibit the central concentra- 
tion of the field (Lubow, Papaloizou, & Pringle 1994). Alterna- 
tive explanations of the filaments have been explored by a num- 
ber of authors (Chudnovsky et al. 1986, Benford 1988, Hey- 
vaerts, Norman, & Pudritz 1988, Lesch & Reich 1992, Rosso & 
PeUetier 1993, Shore & LaRosa 1999). If a pervasive mG field 
is present, it would play a dominant role in the dynamics of the 
interstellar medium near the Galactic center, and would sug- 
gest that the Galaxy was born with a magnetic field > 10^^ G 
(Chandran et al 2000 — this point is discussed further in sec- 
tion 8). 

One of the most pressing questions regarding the pervasive- 
mG-field hypothesis is whether such a strong magnetic field 
could be stably confined. The volume-dominant phase of the 
central region is hot plasma at a temperature T ^ 10^ K which 
fills an elliptical region 150 x 270 pc (FWHM) with a major 
axis that is tilted 20° with respect to the Galactic plane (Ya- 
mauchi et al. 1990). The density of this plasma was estimated to 
be 0.03 - 0.06 cm^^ by Yamauchi et al. (1990) and - 0.3 - 0.4 
by Koyama et al. (1996). The mass-dominant phase of the cen- 
tral region is molecular gas with n > lO''^ cm^^, T ^ 70 K, a 



filling factor in excess of 0.1, and a vertical thickness of ^ 30 
pc (Bally et al. 1988, Morris & Serabyn 1996). If the mag- 
netic field strength is 1 mG, then the magnetic pressure B^ /%% 
is ~ 4 X 10^^ dyne/cm^, which far exceeds the thermal pressure 
of both the hot plasma, /?hot "^4x lO"'" — 4x 10"^ dyne/cm^, 
and the molecular material, /?cioud ^ 10^ "^ dyne/cm^. The ther- 
mal pressure of the ambient medium is thus unable to confine 
mG fields, whether they are ubiquitous or in the form of isolated 
narrow flux tubes. 

The only force capable of confining a pervasive vertical mG 
field is the weight of the molecular material. Interestingly, the 
polarization of dust emission indicates that the molecular ma- 
terial is threaded by horizontal (i.e., || to the Galactic plane) 
magnetic fields (Morris & Serabyn 1996). Because this hori- 
zontal field is perpendicular to the vertical magnetic field in the 
hot plasma, it is difficult for the two types of material to inter- 
penetrate, which suggests that the weight of the molecular ma- 
terial may be able to stably prevent a strong vertical field from 
expanding. Although the molecular clouds are clumpy, it is ex- 
pected that the horizontal magnetic field completely encloses 
the central region. In this paper, the equilibrium of figure 1 is 
proposed as a starting point for investigating the confinement 
of vertical fields at the Galactic center 

To make a compact analytic treatment possible, the stabil- 
ity of this equilibrium is explored under a number of additional 
simplifying approximations. It is assumed that the fractional 
ionization within the molecular clouds is sufficiently large and 
the frequencies of any instabilities sufficiently small compared 
to the neutral-ion collision frequency that ambipolar diffusion 
can be ignored. The partially ionized clouds then behave as 
a single-species fluid to which the magnetic field is frozen, 
or, in other words, as a dense plasma described by the equa- 
tions of magnetohydrodynamics (MHD). The equilibrium is 
taken to be infinite and invariant along the spin axis — the z 
axis in cylindrical coordinates (r, 9,z) — as depicted in figure 2. 
A constant angular velocity £1 is assumed, which removes 
the magneto-rotational instability from the problem (Balbus & 
Hawley 1998). An axial magnetic field and low-density plasma 
are assumed for r < 150 pc, and an azimuthal magnetic field and 
dense plasma for r > 150 pc. The equilibrium density po, pres- 
sure p(j, and field strength Bo are assumed to be proportional to 
powers of r within the dense plasma, and p = 87t/?o/Bo i^ taken 
to be constant and small, which implies that the sound speed a 
and Alfven speed va are proportional to r. The gravitational ac- 
celeration g* is taken to be proportional to r and to be directed 
in the — r direction. If ^* arises from stars, this amounts to as- 
suming a constant mass density of stars throughout the region 
occupied by the plasma. It is also assumed that the plasma is 
supported primarily by rotation, so that va <C Vrot in the dense 
plasma, where Vrot = ^r. Although some of these assumptions 
are extreme, a number of the key ingredients are retained: com- 
pressibility, the Coriolis force, the suspension of dense plasma 
"above" low-density plasma, and the orthogonal orientations of 
the magnetic fields in the low- and high-density plasmas. 

The analysis then focuses on modes which perturb the sharp 
boundary between the high- and low-density plasmas. Such 
perturbations are to some extent stabilized by magnetic tension, 
since short-wavelength disturbances involve significant field- 
line bending of either the axial or azimuthal magnetic fields. 
Uniform rotation also plays a stabilizing role (e.g., Talwar 1960, 
Chandrasekhar 1961, Gilman 1970, Acheson & Gibbons 1978). 
Normal modes of the form ui (^^gl'^^z^+'^^-io?) ^j^j^ ^^ = Q and 



m 7^ are studied analtyically. Such modes avoid bending or 
compressing the axial magnetic field at r < 150 pc and also 
allow compressions of the dense plasma along the azimuthal 
magnetic field that enhance the destabilizing effects of gravity. 
It is shown in sections 2 through 4 that the necessary and suffi- 
cient condition for stability to k^ = Q modes is 

\8\<2\a\a (1) 

when a ^ VA <C iir, where g — g* — ii^r is the effective gravity, 
and a and va are evaluated in the dense plasma. A physical ex- 
planation of equation (1) is given in section 5. Since the effec- 
tive gravity is proportional to the partial support of the plasma 
by the magnetic field, equation (1) can be re- written as an up- 
per limit on the strength of the vertical field [equation (94)]. 
For Galactic-center parameters, the upper limit is ~ I mG. If 
equation (1) is the approximate criterion for the stability of the 
boundary between the high and low-density plasmas in systems 
with differential rotation, a finite extent along the z axis, and 
more general radial profiles — say by evaluating the terms in 
equation (1) near the inner edge of the dense plasma — then a 
pervasive mG vertical magnetic field can be confined at the 
Galactic center. Equation (1) would then also be relevant to 
other accreting systems with strong poloidal magnetic fields. 

Global instabilities with k^ ^ 0, local quasi-interchange 
modes, and instabilities of the two-dimensional plasma of fig- 
ure 1 are discussed in section 6. In section 7, observational 
constraints on the field-strength in the Galactic-center filaments 
are reviewed. The implications of a pervasive mG field at the 
Galactic center for the origin of the Galactic magnetic field are 
discussed in section 8. 

2. DIFFERENTIAL EQUATION FOR k-, = NORMAL MODES IN 
THE DENSE PLASMA 

The plasma is taken to obey the equations of ideal MHD and 
an adiabatic equation of state. 
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where p, p, u, and B are the density, pressure, velocity, and 
magnetic field, g* is the gravitational acceleration, (r, 9,z) are 
cylindrical coordinates, and y is the adiabatic index. An equilib- 
rium is assumed in which all quantities are functions of r alone. 
Equilibrium quantities are denoted with the subscript 0; 

Po^PoW, (6) 

PQ^po{r), (7) 

uo = iirO, and (8) 

Bo = Bo we. (9) 

A uniform angular velocity is assumed, 

Q. = constant, (10) 

and the effective gravity g is defined by the equation 

(11) 



Q}r. 



By virtue of equation (3), 
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and where 
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The sound and Alfven speeds are defined, respectively, by the 
equations 



and 



A- 
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" Po' 

_Bl_ 
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(13) 



(14) 



Small-amplitude disturbances to the equilibrium, denoted by 
the subscript 1, can be decomposed into normal modes of the 
form 

ui=ui(r)e™*'+''^^'-''", (15) 

with analogous equations for Bi, pi, and pi. The equations 
governing these disturbances are obtained in the standard way, 
by substituting u = uo + ui, B = Bq +Bi, etc into equations (2) 
through (5) and discarding all terms containing the product of 
more than one of the perturbed quantities. The analysis is con- 
fined to the case 

kz^O. (16) 

The frequency CO of the disturbance in the frame rotating with 
the plasma is given by 

CO = a- mil. (17) 

From equations (2) and (5), the r and 9 components of equa- 
tion (4), and the 9 component of equation (3), one obtains 
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where 


mf fl2x 

h = 2a&~- 8--), 

r \ r J 
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and 


2 m2fl2 
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Upon substituting equations (18) through (22) into the f com- 
ponent of equation (3), one obtains 
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dr 
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Ds J dr 
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r dr 
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(26) 
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Da = co2 ^. 



(27) 



For isothermal equilibria and Y= 1, equation (25) reduces to the 
k^ ~ limit of equation (2. 13) of Acheson & Gibbons (1978). 
Because k^ = 0, mi-, and Bi^ decouple from the other vari- 
ables and exhibit purely stable behavior: the z component of 
equation (4) gives 

Bu = -^^^^^, (28) 



cor 



and the z component of equation (3) then gives 
Dauiz^O- 



(29) 



3. DIFFERENTIAL EQUATION FOR k^ == NORMAL MODES IN 
THE CONSTANT-p EQUILIBRIUM 

It is now assumed that the equilibrium has constant P, where 

Snpo 



P 



'^O 



(30) 



and that the equilibrium quantities in the dense plasma are pro- 
portional to powers of r: 

Bo{r)=Bx-% (31) 



PoW^px-"?, 



where 



ri 



(32) 
(33) 

(34) 



and ri, B, and p, are constants. It is also assumed that 

g = gx, (35) 

where ^ is a constant. Substituting equations (31) through (35) 
into equation (12), one finds that 

q = 2s + 2, (36) 

VA '^ r, (37) 

fl oc r, (38) 
and 

^ = ^, (39) 

r r'^ 

where 

a=[(<?-2)(l + p)-2]/2. (40) 

It is assumed that a and q are of order unity. 

Using equations (31) through (40) in equation (26), one finds 
that 



F = pHx-'i, 



(41) 



(42) 



where 

^2 , X fl^ / nih \ vl — fl2 

^-^A-^+(2-'?);2(i+75:)+^;^ 

Since Da, D^, h, and H are constant within the dense plasma, 
equation (25) has simple power-law solutions, 

Ulr = X^, (43) 

(44) 
Y = H['-^ + =.^ 1 . (45) 



where 
and 



b'^ + b{2-q)+T = 0, 



9 9 9 \ — 1 



r^D, 



4. DISPERSION RELATION FOR ^z = NORMAL MODES IN 
THE CONSTANT-P EQUILIBRIUM 

The dense plasma is assumed to have an inner edge at r = ri, 
inside of which there is perfectly conducting plasma of vanish- 
ingly small density: 

p^O {forr<n}. (46) 

Three different cases are considered for the behavior at large r. 
(1) the dense plasma has an outer edge at r = r2, outside of 
which is non-magnetized plasma of vanishingly small den- 
sity, (2) there is an outer edge at r = r2, outside of which is 
low-density plasma with an axial magnetic field, and (3) the 
constant- p equilibrium profile is taken to extend to infinitely 
large r. All three cases lead to the essentially the same disper- 
sion relation and stability criterion. In this section, cases (1) 
and (2) are treated; case (3) is discussed in appendix A. 

For uniform rotation, the displacement ^ of a fluid element 
relative to where it would have been in the equilibrium flow is 
given by 

^ = ^. (47) 

CO 

The boundary conditions at each of the boundaries between the 
high and low density plasmas are 

l|n-^||=0, (48) 

and 



and 



(49) 



= 0, 



where n is the unit normal to the perturbed boundary, and 
denotes the jump in / across the perturbed boundary. The La- 
grangian perturbation of a quantity /, denoted 8/, is related to 
the Eulerian perturbation, f\, by the equation 5/ = /i +^ • V/q. 
Equation (49) indicates that the Lagrangian total-pressure per- 
turbation (magnetic plus thermal) is constant across the sharp 
boundary. 
Since 

p-^0 {forr < n orr > r2}, (50) 

a and/or va must be arbitrarily large in the low-density plasma 
in order to maintain pressure balance in the equilibrium state. 
It is assumed that (O is fixed by the dynamics inside the dense 
plasma, so that high-frequency Alfven and sound waves in the 
low-density plasma are not treated. For the comparatively low- 
frequency modes that are considered, the low-density plasma 
instantly rearranges itself to avoid the local compressions 
that would lead to thermal or magnetic-pressure perturbations. 
Within the low-density plasma, the Lagrangian total-pressure 
perturbation thus vanishes, as shown in appendix B. Equa- 
tion (49) then implies that the Lagrangian total-pressure per- 
turbation in the dense plasma also vanishes at the sharp bound- 
aries: 

/,i+^.V/,o+^+^-V^=0 {forr-norr-r2}. 
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(51) 
Using equations (19) and (22), one can rewrite equation (51) 



as 



du\r\ ( ¥2 

A\u\r+A2\x =0 <forx=lorx= — 

ax J [^ r\ 



where 



Ai = -j{Ds + mh), 



(52) 



(53) 



A2 = 



(iS^a^ 



D,vl 



(54) 



The general solution to equation (25) for the constant-p equi- 
librium is a linear combination of the solutions given in equa- 
tion (43), 



Ml, =C1X^1+C2X''2, 




(55) 


where 






^-(f-')W(f-'r-- 


(56) 


and 






'-(f-')V(f-')^-- 


(57) 


Applying equation (52) at x = 1 and at 






_ r2 




(58) 


yields the matrix equation 






/ Ai+A2bi Ai+A2b2 \ f ci \ 


( 


\ 


[ Ai+A2bi x''^^-'''[Ai+A2b2] )\C2 ) 


"lo 


(59) 


Non-trivial solutions require that 






(Ai +A2fel)(Ai +A2fe2)(x2'"^' - 1) = 


= 0. 


(60) 


One type of normal mode results from setting 


^2 ~ 


l,or. 


equivalently. 






b2-bi = ^ /"' ,, n = ±l,±2, 




(61) 



\n{r2/ri] 
Equations (56), (57) and (61) yield the dispersion relation 

T-pn = Q, (62) 

where 

%n 



Pn 



n2 



ln(r2/ri) 



(!-')■ 



(63) 



The Lagrangian total-pressure perturbation of these modes van- 
ishes iiir — r\,r — r2, and at « — 1 discrete values of r between 
ri and r2- Note that n = in equation (61) does not correspond 
to an eigenmode, since it implies by = b2, cy = — C2 from equa- 
tion (59), and thus u\r = Q. 

A second type of normal mode results from having either 

Ai+A2foi==0 and (64) 

C2 = 0, (65) 



or 



Ai+A2fe2==0 and (66) 

ci == 0. (67) 

In either case, one sets b ~ — A1/A2 and substitutes into equa- 
tion (44) to obtain 

^4 



0=^(Z)s 



nhy 



^iDs 



fl^CO^ 



DaDs - h^ 



D,(vi-fl2) 



(68) 



The Lagrangian total-pressure perturbation of these modes van- 
ishes everywhere — they are the gravest modes of the system 
for any given m. The right-hand side of equation (68) is pro- 
portional to Ds, but the Ds = roots are spurious; in deriving 
equation (25) it was assumed that D^ ^ 0, and it can be verified 



that Ds = does not lead to normal-mode solutions satisfying 
the boundary conditions. 

Equations (62) and (68) can be rewritten in the form 

= D^Dr + S, 



where 



CO 



Dr = (0 - (0 



4mgD, 



m 
~G 



Lf^v 



2„2 
A 



r4 



Cv? 



and 



Pn 



m^ - 1 
-1 



for the roots of equation (62) 
for the roots of equation (68) 



(69) 
(70) 
(71) 

(72) 
(73) 



The value of S is different for equation (62) and equation (68). 
For equation (62) 

4mgD.(ji^ m^K?- f g^ C^^^a 



S = 



rG 



Co2(^ + Q+m^(2(on-^) 



(74) 



whereas for equation (68) 
S 



2 2 



-2m(2coi2 



CO 



1-- 



2 2 

ma 



m^{m^ — l)fl^ 



,_'ng_\ 



4-mgQ.(i)^ m^(ii^{AQ.^a^ +g 



(75) 



rG Gr^ 

It is now assumed that the dense plasma satisfies the low-p 
condition 

fl^VA, (76) 

and is supported against gravity primarily by rotation, which 
implies that 

VA < Vrot, (77) 



where 



Vrot= \^r\. 



(78) 



Equations (76) and (77) imply that S is negligible compared to 
the dominant terms in DgDu for all values of co, m, and p„. The 
approximate roots of equation (69) are then obtained by setting 
DeDR = 0. When D^ = 0, 



CO: 



±j4ii2 + ^vyr2. 



When Dr = 0, 



CO: 



2mgQ.r^ ± mvAd^i'^^^a^r^ - g^r^ + ^a\\) 



AD?r^ 



-Wk 



(79) 



(80) 



The errors in equations (79) and (80) vanish as (a/va) -^ and 
(vA/vrot) -^ 0. Comparisons between equation (80) and numer- 
ical solutions of equations (62) and (68) are given in figure 3. 

When fl/vA and VA/vrot are small (just how small will be 
seen below in the numerical solutions), both the approximate 
roots given in equations (79) and (80) and the exact roots of 
equations (62) and (68) are real for all m and n if and only if 

\g\<2\a\a. (81) 

[To be precise, the exact stability criterion is slightly differ- 
ent from equation (81) due to small corrections associated with 



nonzero a/vA and vp^/vmv] It can be seen directly from equa- 
tion (80) that the approximate roots are real for all m and n if 
and only if equation (81) is satisfied; that the exact roots are 
also real can be seen from the following argument. Let 

D = D^Dr + S. (82) 

The four roots of the exact dispersion relation D — Q are real 
when the plot of Z)(co) passes through zero four times, as in 
figure 4. When a <C va and va ^ Vrot, D always has two 
local minima at co ~ ±\jG/2 at which D ~ —0^/4. Also, 
D -^ +00 as |co| ^ +00. The question of stability turns upon 
whether D is non-negative for some value(s) of CO in the in- 
terval (— ■\/G/2, yjG/2). When a ^ va and va <C Vrot and 
equation (81) is satisfied, the local maximum of 0^.0^ at 
<M = ©max — 2mgQ./rG G (— ■\/G/2, yjG/2) is positive. More- 
over, 5(c0max) is negligible compared to the value of D^D^ at 
Wmax (unless the plasma is extremely close to marginal sta- 
bility). Hence, the value of D is also positive at co^ax G 
(— ■\/G/2, yjG/2), and the four exact roots are real. On the 
other hand, when equation (81) is not satisfied, the local max- 
imum of DeDji at CO = COmax is negative, D^Dn is negative for 
all CO in the interval {-yjWf2, ./gJ2), \S\ < \D^D^\ for all CO 
in the interval (— ■\/G/2, yjG/2), D^D^ + 5 is negative for all 
CO in the interval (— ■\/G/2, ^G/2), and two of the exact roots 
have nonzero imaginary parts. 

The \m\ — 1 roots of equation (68) are a special case, since 
equation (80) gives a double root at CO = g/2Q.r. In this case, the 
local maximum of D^D^ occurs at CO = g/2Q.r, and DgDji_ = 
at the local maximum. It is thus not clear whether the local 
maximum of D^Dj^ + S has a positive or negative value. In ap- 
pendix C, a new approximate factorization of equation (68) is 
given for |m| = 1 and it is shown that the |m| == 1 roots are 
always stable when a <C va and va <C v^ot, whether or not equa- 
tion (81) is satisfied. 

Figure 5 summarizes a stability analysis based upon numeri- 
cal solutions of equations (62) and (68) for values of a/va and 
VA/vrot ranging from 0.01 to 1. A combination of parameters 
(o/va) and (vA/vrot) is said to be stable if the imaginary part 
of each of the numerically determined roots is less than lO^'*' 
times the real part of the same root for all values of m from 1 
to 100, and in the case of the roots to equation (62), for all val- 
ues of n from 1 to 500, where ln(r2/ri) is taken to be 5, and q 
is taken to be 5. The roots are determined using an eigenvalue 
method (Press et al 1992) at double precision. Equation (81) 
describes the stability of the exact roots quite accurately pro- 
vided VA/vrot < 0.5. For larger values of va/vtoi, equation (81) 
underestimates stability. 

For any given values of a, Q., and g, the most unstable mode 
has zero Lagrangian total-pressure perturbation everywhere and 
^ = m^ — 1. Upon setting d03i/dm — 0, where CO/ is the imagi- 
nary part of CO, one can show using equation (80) that the max- 
imum growth rate Ymax occurs at 
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(84) 



where / is a function that increases monotonically from to 
1 as its argument Q?a^/g-^ is decreased from 0.25 (marginal 
stability) towards 0. Because of equation (39), Ymax ~ va/^ 



when the plasma is not very near marginal stability. When and using equation (18) to eliminate Bir, one obtains 



\g\ > 2\D.\a, modes with ^ — m — 1 as well as modes with 
^ — Pn + rn^ — ^ are unstable only when C, < ^crit, where 



r^{g^~4-Q.V) 



2„2 



a^v 



(85) 



5. DISCUSSION OF STABILITY CRITERION AND 
EIGENMODES 

In the absence of background flow, the stability of an equi- 
librium to a small-amplitude displacement ^(x) is determined 
by whether ^(x) leads to an increase or decrease in 5W, where 
8W is the sum of the perturbed magnetic, internal, and gravi- 
tational potential energies of the plasma (Bernstein et al 1958). 
If dW is positive, then the mode is stable — if the displacement 
increases with time initially, the initial bulk-motion kinetic en- 
ergy present in the mode is converted into these other types 
of energy, and the plasma proceeds to oscillate stably. If dW 
is negative, then the magnetic, internal, and gravitational ener- 
gies of the plasma can be tapped to drive an instability. In an 
equilibrium with flow, the kinetic energy of the flow is altered 
by plasma displacements, and stability is no longer determined 
solely by the increments to the magnetic, internal, and gravita- 
tional potential energies. The perturbed kinetic energy depends 
upon both the displacement ^(x) and its time derivative dt,/dt, 
which gives rise to the phenomenon of overstability, in which 
both the real and imaginary parts of (O are nonzero, and which 
makes the stability analysis more complicated than in the static 
case (Frieman & Rotenberg 1960). 

In this section, the physical basis of the stability criterion is 
explained using the properties of the eigenmodes. For any set 
of plasma parameters, the most unstable k^ = modes are those 
described by equation (68) [as opposed to those described by 
equation (62)], for which the Lagrangian total -pressure pertur- 
bation vanishes everywhere. It is these modes that will be dis- 
cussed in this section. The dispersion relation given by equa- 
tion (68) has two branches, described approximately by equa- 
tions (79) and (80). Equation (79) is an epicyclic-Alfvenic 
branch. When m ~ (1), one finds that co = ±2ii, uiq ~ T'^ir, 
and b ~ — (±ot+ l)a^/v\. These circularly polarized epicyclic 
modes, which involve primarily a balance between inertia and 
the Coriolis force, are stable for all plasma parameters: when 
a fluid element in the dense plasma of figure 2 tries to sink to- 
wards r = 0, the Coriolis force pushes it in the +9 direction (for 
positive i2); as it moves in the +9 direction, the Coriolis force 
pushes it back towards larger r. For rrp- 3> Vj(,(/v^, the modes 
are stable Alfven waves, with co = ztrnvA/r, \uiq\ ^ \uir\, and 
b ~ — fl^/v^. The dominant force across field lines for these 
modes is magnetic tension. 

Equation (80) is the Rayleigh-Taylor/sound-wave branch of 
the dispersion relation. It is shown in appendix D that the eigen- 
modes corresponding to equation (80) satisfy JMiel ^ |Mir|, 
which indicates that these modes involve comparatively little 
field-line bending. Form^ o(l), a)<Cii. In the frame rotating 
with the plasma, the momentum equation at first order in pertur- 
bation amplitude can be written — /copoui = — 2poiiz x ui + . . ., 
where only the Coriolis force Fc has been written explicitly 
on the right-hand side. Since co <C ii and Ui is not along z, 
some force must approximately cancel the Coriolis force. For 
m ~ ( 1 ), the near cancellation is provided by the gravitational 
force F^. Upon defining F^, || — Fg ■ B/Bq and Fc,ll — Fc • B/Bo, 



pomguir 
F„\\ = , and (86) 

Fcjl = -Ipo^uir. (87) 

It is shown in appendix D that pi ~ pomMie/cor (that is, the 
density perturbation is dominated by the compression associ- 
ated with Uiq), so that 

PomguiQ 



while 



^ cor 

r-Fc =2poiiMie. 



(88) 

(89) 

The phase velocity of a low-m mode in the 9 direction, Vphase = 
COr/m, is set by the requirement that Fc and F^ approximately 
cancel, which yields Vphase — §/2ii, or CO ~ mg/2D.r. This 
phase velocity is in the eastwards direction when g > — i.e., 
in the direction of rotation (co/mi2 > 0). Because of the addi- 
tional forces in the plasma, which are much smaller than either 
Fg or Fc, the exact phase velocity will be either a little larger 
than |^/2i2|, in which case \Fg\ < \Fc\, or a little smaller than 
|^/2ii|, in which case \Fg\ > \Fc\- 

It will prove useful to consider the radial momentum equa- 
tion in somewhat greater detail. It is shown in appendix D that 
\b\ <C 1 for the eigenmodes of equation (80) with m ^ 0(I). 
This is related to the fact that the Lagrangian perturbation of 
the magnetic pressure, which is <=<: duir/dr, is cancelled by the 
Lagrangian perturbation of the much smaller thermal pressure. 
(Although |Mie/Mir| is large, it is not large enough to enable \b\ 
to be > 1.) With the use of the relation \b\ <C 1, the vanishing 
of the Lagrangian total-pressure perturbation, equations (18) 
and (19), and the equation pi = — poV • ^ — ^ • Vpo (Frieman 
& Rotenberg 1960), the radial component of the momentum 
equation can be written 

mg \ iuirPojm^ - l)v\ 



(COpoMlr 



PoMie 



2i2co 



(90) 



CO \ r J cor'' 

The (poMie/co)(2i2co — mg/r) term is the radial component of 
what will be called the net gravitational force, F^, +Fc. The 
—iuirPo{m^ — 1)i'a/'^''^ '^srm is the effective magnetic tension, 
which acts to restore a radially displaced fluid element in the 
dense plasma to its initial position. In order for a mode with 
CO ~ mg/2Q.r to propagate, the radial components of the effec- 
tive tension and the net gravitational force must approximately 
cancel; otherwise, solving for CO yields |co| > VA/r, which is 
^ mg/2Q.rfoi m ^ 0(1). 

Because gravity and the Coriolis force nearly cancel, the low- 
m eigenmodes of equation (80) do not enjoy the full stabiliz- 
ing effect that the Coriolis force has on epicyclic modes, nor 
the full destabilizing effect that gravity has on Rayleigh-Taylor 
modes in a static plasma. The physics behind the stability cri- 
terion can be understood by considering a mode with m ^ O ( 1 ) 
and |vphase| slightly less than |g/2ii|, which implies that \Fg\ 
is slightly bigger than \Fc\ and that the net gravitational force 
acts in the same direction as gravity. If |vphase| < o, then the 
dense plasma reacts in a way that is similar to how it would 
react if Vphase were 0: the maximum density enhancements oc- 
cur at those points towards which the parallel component of the 
net gravitational force converges — the troughs of the perturbed 
magnetic field lines, as in figure 6. The induced pressure force 
opposes the parallel component of the net gravitational force. In 
the radial direction, the net gravitational force and the effective 
magnetic tension oppose each other and can therefore approx- 
imately cancel, as required for a mode with CO ~ mg/2Q.r. On 



the other hand, if |vphase| > a, then sound waves can not prop- 
agate information about the net gravity upstream of the wave. 
At any instant in time, fluid elements at the troughs of the per- 
turbed field lines can not tell that the net gravitational force is 
trying to compress them, and the plasma reacts inertially to the 
parallel component of the net gravitational force. Plasma is ac- 
celerated towards positive 9 at those points where Bir/Bo < 0, 
and the resulting azimuthal compression causes the density en- 
hancements to occur at the crests of the field lines when (O is 
real, as in figure 6. [This point can also be seen with the aid 
of equation (D18) in appendix D.] In this case, the net gravi- 
tational force in the radial direction acts in the same direction 
as the effective magnetic tension, providing much more radial 
force than is consistent with co ~ mg/2Q.r, and no stable mode 
at approximately that frequency exists. 

Overstable modes, however, do exist. One way of thinking 
about this is the following. Suppose a low-m mode were to start 
propagating with a phase velocity greater than a and peak den- 
sity enhancements at the crests of waves (as would be the case 
for real co). The radial component of the net gravitational force 
and magnetic tension would then provide much more force than 
required to return the displaced plasma to its initial position. A 
fluid element displaced a distance t,ro towards +r would in re- 
sponse pass through its initial position and then move a distance 
> t,rQ towards — r, leading to an oscillation of growing ampli- 
tude. If (0/ 7^ 0, where (0/ and cor are the imaginary and real 
parts of CO, respectively, then the phase relation between p i and 
uii- is not as depicted in figure 6. When m ^ 0(1), |cOff| 3> |cO/| 
as in equation (80). Fg\\ and Fqh are then effectively sine waves 
of equal amplitude that are nearly 180° out of phase, and their 
sum is about 90 degrees out of phase with each of them. The 
locations at which pi obtains its peak magnitude are then ap- 
proximately the steepest points on the perturbed magnetic field 
lines, in phase with Bir and ui,-. The peak magnitude of the 
radial component of the net gravitational force is then about 90 
degrees out of phase with p i and coincides with the troughs and 
crests of field lines. The radial component of the net gravity is 
then approximately cancelled by the effective magnetic tension. 

When \g\ > 2|i2|fl, modes are stabilized when ^ is increased 
above ^crit because pressure begins to dominate over gravity 
in the dynamics along the magnetic field lines, changing the 
modes from Rayleigh-Taylor instabilities into something akin 
to sound waves. It is shown in appendix D that \Fg\\ \ > \Fp\\ \ 
when ^ < ^crit and 1^1 > 2|i2|fl, where Fpn is the component of 
the pressure force along the perturbed magnetic field lines in 
the dense plasma. It is also shown that as ^ is increased towards 
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l.ForC>Ccrit, l^p. 



> l^..| 



for all m with 



one exception: ^crit < C < ^i''^ l'^^\ on the branch of the dis- 
persion relation with ±i2 < 0. This exception corresponds to 
the segment of the plot of cor in the upper-left panel of figure 3 
just to the right of ^crit where the curve bifurcates, between the 
bifurcation point and where the plot of cor intersects cok = 0; the 
plot of (Or bifurcates at m — 49, and the exception corresponds 
to the smaller value of co at m := 50. 

To obtain further insight into the stability criterion, it is help- 
ful to consider some limiting cases. Of course, if there is neither 
rotation nor magnetic field, the equilibrium of figure 2 is unsta- 
ble to Rayleigh-Taylor modes as in the classic calculation of a 
dense fluid suspended above low-density fluid in slab geometry 
(Rayleigh 1883). If the orthogonal magnetic fields of figure 2 
are present but the plasma isn't rotating, the situation is simi- 
lar to the stationary plane-parallel problem solved by Gratton, 



Gratton, & Gonzalez (1988) and Gonzalez & Gratton (1990). 
These authors discussed a number of cases, including that of a 
plasma at >■ > suspended above vacuum at y < with grav- 
ity pointing in the —y direction, all equilibrium quantities being 
functions of y alone, the magnetic field in the plasma pointing 
in the x direction, and the magnetic field in the vacuum point- 
ing in the z direction. They found that the orthogonal magnetic 
fields act to stabilize the plasma, since the interface at y = 
can not be perturbed without bending either the magnetic field 
in the plasma or the magnetic field in the vacuum. Compress- 
ibility, however, acts to destabilize the plasma; if the plasma is 
strongly compressed along the magnetic field, the gravitational 
force acting on regions of enhanced density can overcome mag- 
netic tension. When the wave vector k of the perturbation in the 
X — z plane is oriented along the magnetic field in the plasma at 
y > 0, magnetic tension on its own is not sufficient to stabilize 
small-scale modes: as a ^ 0, the maximum unstable wavenum- 
ber approaches o°. Since thermal pressure limits the compress- 
ibility of the plasma along the magnetic field, magnetic tension 
and pressure together stabilize modes with k^ > fe^^^j, where 
^crit ^ ^^(^A + fl^)/fl^v'^- It can be shown that the condition 
1^1 > ^crit corresponds to pressure dominating over gravity in 
the dynamics along the perturbed magnetic field lines, in anal- 
ogy to the stabilization of the modes described by equation (80) 
when ^ > ^crit- On the other hand, sufficiently large-wavelength 
modes (|A:| < k^xxd ^e always unstable in the magnetized slab, 
whereas all k^ — ^ modes are stable in the rotating magnetized 
cylinder when \g\ < 2\Q.\a. 

If there is rotation but no magnetic field, the plasma has a set 
of normal modes analogous to those described by equation (68) 
(i.e., k^ — and zero Lagrangian pressure perturbation every- 
where), with frequencies 
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These modes are circularly polarized (uiq = ±Mir), stable at 
long wavelengths (|m| < Q?r/g), unstable when m > Q.^r/g 
in equation (91), and unstable when m < —Q?'r/g in equa- 
tion (92). As \m\ -^ °o, the growth time of unstable modes 
(^ Xj \fmgjr) is much shorter than the rotation period, and 
rotation can be neglected. The introduction of an azimuthal 
magnetic field into the dense plasma has a stabilizing influence 
insofar as magnetic tension eliminates these large- |m| instabili- 
ties with JMiel = \u\r\. 

The magnetic field in the equilibrium of figure 2, however, 
can also be destabilizing: as the field strength is increased, the 
plasma becomes increasingly supported by magnetic pressure, 
thereby increasing the effective gravity g. With the use of equa- 
tion (39), the stability criterion can be rewritten as an upper 
limit on the Alfven speed: 

^ 'flVrot (93) 



A< 



a 



(where a, va, and Vrot are evaluated in the dense plasma). In 
the low-p limit, pressure balance requires that the magnetic- 
field strength Byert in the low-density plasma at r < ri be ap- 
proximately equal to the field strength in the dense plasma near 
r ~ ri, B. Assuming a dense plasma of pure hydrogen with 
q = 5,J^ 5/3, p < 1, Vrot = 150 km/s [from the Galactic ro- 
tation curve at r = 150 pc (Ruzmaikin & Shukurov 1981)], and 



r = 70 K and n~ 10"^ cm ^ near r = ri, one can rewrite equa- 
tion (93) as 

Bvert<l.lmG. (94) 

Although it is tempting to infer from equation (94) that mG 
vertical magnetic fields at the Galactic center can be stably con- 
fined by a ring of horizontal magnetic fields threading molecu- 
lar material, further study is needed to extend the present anal- 
ysis to more realistic equilibria. 

6. OTHER TYPES OF INSTABILITIES 

Sections 2 through 5 treated the global stability of the equi- 
librium pictured in figure 2 to k^ = modes. This section ad- 
dresses three other classes of instabilities: (1) local interchange 
and non-axisymmetric quasi-interchange modes in a uniformly 
rotating plasma with azimuthal magnetic field, (2) global in- 
stabilities in the ID equilibrium of figure 2 with k^ ^ 0, and 
(3) quasi-interchange and magneto-rotational instabilities in a 
two-dimensional differentially rotating equilibrium such as the 
one pictured in figure 1 . 

Local modes in a plasma with an azimuthal magnetic field 
that rotates uniformly about the z axis have been treated by 
a number of authors. For the constant-p equilibrium pro- 
file of section 3, local axisymmetric interchange instabilities 
are suppressed primarily because of the strongly stabilizing 
specific-angular-momentum profile when a ^ va ^ ^rot (New- 
comb 1962, Acheson & Gibbons 1978, Rogers & Sonnerup 
1986, Ferriere et al 1999). Such modes are further stabilized 
since the specific entropy and magnetic flux per unit mass of 
an azimuthal flux loop B/pr are taken to increase with r in 
the constant-p equilibrium of section 4. The results of Fer- 
riere et al (1999) on a wide class of one-dimensional equilib- 
ria can be used to show that local non-axisymmetric quasi- 
interchange modes in a uniformly rotating plasma with az- 
imuthal magnetic field are stable when two conditions are met 
assuming a <IC va ^ Vrot: (1) the "magnetic Brunt-Vaisala" fre- 
quency COo must be real, and (2) |^| < 2|ii|fl (to lowest order 
in a/vA and VA/vrot). For the constant-p equilibrium of sec- 
tion 4, ©0 — qg/'^r and the first condition is met. The sec- 
ond condition is the same as equation (81), and thus when the 
k^ — Q modes are stable, so are the quasi-interchange modes. 
The quasi-interchange modes have co ~ mg/2Q.r, and are stable 
(provided coo is real) when their phase velocity along the mag- 
netic field is less than the sound speed, just as for k^ = modes 
with m^ 0{1). 

Global modes in the equilibrium of figure 2 with k^^O and 
m ^ ( 1 ) are analogous to the external kink instability in static 
plasmas. An m—\ mode is illustrated in figure 7 [figure 9.23 
of Freidberg (1987)]. Because the strength of the azimuthal 
field increases where the field lines are compressed, the Lorentz 
force of the azimuthal field acts to increase the perturbation. On 
the other hand, the magnetic tension of the axial magnetic field 
opposes the perturbation. In the static reversed field pinch, in 
which Be ^ B^, the azimuthal magnetic field wins out over the 
axial field and the mode is unstable when the axial wavelength 
fc-T ' is greater than roughly the radius of the plasma cross sec- 
tion (Friedberg 1987). How this instability mechanism affects 
the stability criterion of the rotating equilibrium of figure 2 is 
left for future work. 

Newcomb (1961) showed that the convective stability crite- 
rion for a compressible stratified magnetized slab when k\\ = 
is less stringent than the stability criterion in the limit kn -^ 0, 



where kn is the component of the wave vector along the mag- 
netic field. This is because when kn — 0, Newcomb's type 2 
modes, which involve motion primarily along the magnetic 
field, become pure translations which are neutrally stable (co = 
0). The fell = stability criterion is thus the one that applies 
to Newcomb's type 1 modes, which involve motion primarily 
across field lines. On the other hand, when kn is small but 
nonzero, the type 2 modes allow destabilizing compressions of 
the plasma along the magnetic field and can be unstable even 
when type 1 modes are stable. The distinction between kn —0 
and fell — > in Newcomb's work is not directly analogous to the 
distinction between fe^ = and fe^ ^ in the equilibrium of fig- 
ure 2, in that fen corresponds to m since it is the dense plasma 
that dominates the dynamics of the unstable modes. 

The two-dimensional equilibrium of figure 1 is subject to in- 
stabilities not present in the one-dimensional infinite-cylinder 
equilibrium. If the vertical (along z) stratification does not sat- 
isfy the Schwarzchild criterion for a static non-magnetized gas, 
then the plasma is unstable to local non-axisymmetric quasi- 
interchange modes with fe^ ^ fe^ and fe^r/m ^ 1 (Terquem & 
Papaloizou 1996). The displacements in such modes are es- 
sentially parallel to the spin axis, and thus the Coriolis force 
vanishes, preventing the type of rotational stabilization that 
occurs for quasi-interchange modes in the one-dimensional 
infinite-cylinder equilibrium. If there is differential rotation 
with d\Q.\/dr < 0, the plasma is subject to the magnetoro- 
tational instability (MRI) (Ogilvie & Pringle 1996, Terquem 
& Papaloizou 1996, Balbus & Hawley 1998). Global three- 
dimensional numerical simulations suggest that turbulence ex- 
cited by the MRI will be present quite generally in astrophysi- 
cal disks (Hawley, Balbus, & Stone 2001). It is thus interesting 
to ask how the stability of the sharp boundary in figure 1 is af- 
fected by differential rotation in a quiescent equilibrium flow as 
well as by differential rotation with turbulence. These questions 
are left for future work. 

7. ESTIMATES OF FIELD STRENGTHS IN RADIO FILAMENTS 

Lang et al. (1999b) found that if the magnetic energy density 
in the Northern Thread is equal to the relativistic-particle en- 
ergy density, and if the relativistic protons and electrons have 
equal energies, then the field strength is ~ 140/jG. Comparable 
equipartition field strengths have been found for the Radio Arc 
(Tsuboi et al. 1986, Yusef-Zadeh & Morris 1987b, Reich 1994). 
However, the equipartition estimate may lead to an underesti- 
mate of the field strength, particularly since the geometry of the 
filaments allows relativistic particles to stream directly out of 
the Galactic-center region. 

All of the radio filaments that have been studied in detail 
near the Galactic center appear to be interacting with at least 
one randomly moving molecular cloud, and yet most are unde- 
formed at the sites of interaction (Morris & Serabyn 1996). If 
the magnetic field were dynamically insignificant, one would 
expect the turbulent motions of the clouds to cause turbulent 
bulk motions in the hot plasma that would deform the frozen- 
in magnetic field lines. The lack of deformation gives a lower 
limit on the field strength in the filaments that depends upon 
how many independently moving clouds (or clumps within a 
cloud) interact with a single filament. 

If a vertical filament interacts with a single molecular cloud 
at a single point along its length, then the lower limit on its field 
strength is fairly small. Since the vertical field Byert threads 
the very tenuous hot plasma, there is very little mass to anchor 
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the ends of the vertical field line above and below the Galac- 
tic disk. Any deformation in the field line in the galactic plane 
by an interaction with a molecular cloud is carried away along 
the field line at the Alfven speed of the hot medium, v^ hot- In 
order for the visible deformations in the filaments to be small, 
the Alfven speed in the hot medium must be much larger than 
the turbulent speed of the the clouds Vdoud with respect to the 
ends of the filament. The angle at which a field line is bent 
will be ^ Vcioud/vA,hot, as depicted in figure 8. The condition 
VA,hot > Vcioud can be written 



1/2 / Vcloud , „ 



(95) 



If 



where nt,ot is the number density of the hot plasma in cm^^ 
Vcioud = 30 km/s, and nhot = 0.3, then Byert 3> 8/jG. 

As pointed out by Yusef-Zadeh & Morris (1987a) and Morris 
& Serabyn (1996), if a filament interacts with clouds at more 
than one point along its length, as happens in the Radio Arc 
(Morris 2001), then the lower limit on the field strength may 
be considerably greater than in equation (95). For certain cloud 
shapes, a filament could conceivably avoid the appearance of 
deformation by rapidly slipping around the moving clouds. For 
some cloud shapes, however, this would not be possible, and 
to avoid deformation a filament would have to dynamically op- 
pose the force applied by the clouds. If this requires that the 
magnetic pressure in the filament exceed the ram pressure in 
the clouds, then Byert > 1 mG (Yusef-Zadeh & Morris 1987a). 
If a single cloud consists of independently moving clumps, then 
a similar argument would apply to a filament interacting with a 
single cloud (Morris & Serabyn 1996). 

8. CONCLUSION 

Because of the simplifying assumptions made in this paper, 
further work is needed to determine the extent to which the sta- 
bility criterion given by equation (81) applies to more realistic 
equilibria, and to extend the analysis to modes with k^ ^ 0. If 
equation (81) does describe, even approximately, the stability of 
the boundary between the high and low-density plasmas in sys- 
tems with differential rotation, a finite extent along the z axis, 
and more general radial profiles — say by evaluating the terms in 



equation (81) near the inner edge of the dense plasma — then it 
appears that a pervasive mG vertical magnetic field can be con- 
fined at the Galactic center Equation (81) would then also be 
relevant to other accreting systems with strong central poloidal 
magnetic fields. 

A pervasive mG vertical magnetic field in the inner 150 pc 
of the Galaxy would have important implications for the ori- 
gin of the Galactic magnetic field. If the vertical field origi- 
nated from radial inflow over the lifetime of the Galaxy, then 
the strength of the average vertical magnetic field in the Galac- 
tic disk at the time of Galaxy formation, B,-, can be estimated by 
magnetic-flux conservation. If the average mass accretion rate 
into the Galactic-center region is 0.3 Mo/yr (Morris & Serabyn 
1996), and if the surface density of gas in the disk has remained 
constant over the lifetime of the Galaxy, then all the gas that 
was within roughly 10 kpc of the Galactic center at the time 
of Galaxy formation has been accreted into the Galactic-center 
region. If the magnetic field was frozen in to the interstellar 
plasma, then the vertical magnetic flux that initially threaded 
the central 10 kpc of the disk, <I>, = B,7i(10 kpc)^, is now con- 
centrated near the Galactic center. If there is a pervasive verti- 
cal magnetic field at the Galactic center today, it is unlikely that 
it contains regions of oppositely directed vertical fields, since 
the parallel geometry of field lines would permit rapid destruc- 
tion of the oppositely directed fields by turbulent diffusion and 
reconnection. The current magnetic flux through the Galactic- 
center region corresponding to such a pervasive field, <!>/, can 
thus be estimated by multiplying field strength times area. If 
there is a mG vertical field throughout the central 150 pc of 
the Galaxy, then ^f = 7t(150 pc)^ mG. Setting 4>, = 4>/^ yields 
Bi ^ 10^^ G (Chandran et al 2000), which would be consistent 
with a turbulent dynamo operating in the protogalaxy (Kulsrud 
et al. 1997). 
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APPENDIX 
DISPERSION RELATION WHEN CONSTANT-p AZIMUTHAL-FIELD EQUILIBRIUM EXTENDS FROM r = ri TO r = oo 

In this case, the boundary conditions are that the Lagrangian total-pressure perturbation vanishes at r = ri and that the behavior as 
r ^ CX3 is such that the kinetic energy of the mode is finite. As in section 4, the boundary conditions can be satisfied in two ways. The 
first is for both c\ and ci to be non-zero. The boundary condition at r = ri is then given by the first row of equation (59), and a finite 
kinetic energy requires that the radicals in equations (56) and (57) be purely imaginary; if they are real or complex, then the r''' term 
causes the kinetic energy to diverge, where JT7. is taken to be the root with a positive real part. The resulting dispersion relation is 

1 )]^ + A^, where A^ is any positive real number The second way is to have 
(0.5^ ^ 1) ^ ^Z (0.5(7— l)^^r^ which leads to equation (68) and corresponds to a mode in which the Lagrangian 
total-pressure perturbation vanishes everywhere. 

DIFFERENTIAL EQUATION FOR fc, = NORMAL MODES IN THE LOW-DENSITY AXIAL-FIELD REGION 

The differential equation for A:^ = normal modes in a uniformly rotating plasma with axial magnetic field [Bq = BQ,(r)z\ is 



the same as equation (62), but with /?„ replaced by [(0.5^ 
u\r °^ r^ with b 
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When g* —0 (i.e., g — —Q?r), equation (Bl) reduces to the constant-ii Hmit of equation (3a) of Spies (1978). In the limit po -^ 0, 
with po constant, equation (Bl) reduces to 
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with solutions 

with 

l = -\±m. 

Writing the Lagrangian pressure perturbation as hp and the Lagrangian perturbation of B^ as hB,, one finds that 
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which vanishes as p ^ for fixed co, where co is set by the dynamics inside the dense plasma. 

APPROXIMATE FACTORIZATION OF DISPERSION RELATION [EQUATION (68)] FOR |m| = 1. 

For the |m| = 1 roots of equation (68), the dispersion relation can be written as 
where Dr = (co - (0+)((0 - (0_), and 



COi 



2Q. 1 r 



1± 1 



4fl2n2 



(CD 



(C2) 



The value of 5 can be recalculated and is different from the values in equations (74) and (75). As in the case \m\ ^ 1, when a <C va 
and Va ^ Vrot, D^D^ + S has local minima at co ~ ±-\/G/2 at which D^D^ + 5 ~ — G^/4, and stability again turns upon whether 
DbSr + 5 becomes non-negative in the interval CO e (— ■\/G/2, ■\/G/2). The value of De^R has a local maximum at co = (co+ + co_)/2 
at which D^Dr is always positive. It can be shown that S <C DqDr at CO = (co+ + co_)/2, and thus D^D^ + 5 is always positive at 
CO — (co+ + C0_)/2, which guarantees that the four exact roots are always real. 

PROPERTIES OF THE EIGENMODES 

This appendix treats the eigenmodes corresponding to equation (68), for which the Lagrangian total-pressure perturbation vanishes 
everywhere. All of the results of this section assume that 

a < VA < Vrot, (Dl) 

and a number of the equations to be presented assume that the plasma is not very near marginal stability, i.e.. 



-4nV 



> — , and 

VA 



^2 + 4n2fl2 

Upon substituting equation (80) into equation (24), one finds that 

^2 



-4n2fl2 



g^+4QP-a^ 



> 



VA 
Vrot 



D, 



^2q2 



la^a^G-g^T 



gVA 



VI 



(D2) 



(D3) 



where the ~ sign has the same meaning as in equation (80), that the error in the equation vanishes as a/va — > and VA/vrot -^ 0. 
When CO is real, (fl2G — g-^) > 0, D^ is real, and D^ < 0. This implies that all stable modes on the Rayleigh-Taylor branch of the 
dispersion relation [equation (80)] travel more slowly than sound waves. Upon substituting equation (80) into equation (23), one 
finds that 



/z~± 



mvA 



vc 



r2G 



2Q.y/a^G-g^T 



gVA 



v~c 



ma 



(D4) 



Given equation (D2), the ma^ /r^ term on the right-hand side is negligible compared to the other term for |m| > 1. It can be shown 

that 

h 



£>s 



>1 



(D5) 



for Iml > 1. 
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When 

1 < m <C Vrot/fl 

and equation (D2) is satisfied, it can be shown that jco^a^l <C li'sV^I, 

a mh 



viD,- 



and, for unstable modes with a G — g <0, 



H-- 



CG 



va\ 8^- 4i22fl2 ■ 



(D6) 



(D7) 



(D8) 



The value of \b\ in equation (D8) is <C 1 when \m\ ^ 0{1) and increases monotonically with increasing ^, reaching a maximum value 
of 



1^1 max — 2 — ^ 



(D9) 



at ^crit, the largest value of ^ for which modes are still unstable [equation (85)]. When m ^ (1), |fe| <C 1 for stable modes as well 
Given equations (D2) and (D6), it can be shown that \{ma'^/r)duir/dr\ <C \huir\ and 

ihuir 



Mie 



Equation (D5) then implies that 

Equations (D7) and (DIO) can be used to show that 



Miel > \uir\- 

imuiQ 



V u 



Since | V • ui | ^ \u\r/r\, the relation (Frieman & Rotenberg 1960) 

pi = -ypoV • ^ - ^ • V/50 
can be used to show that 



;5i~-y/5oV-^: 



ia pomhuir 



Similarly, the relation pi = — poV • ^ — ^ • Vpo (Frieman & Rotenberg 1960) can be used to show that 

pomuiQ ipomhuir 



Pi 



cor 



COrDc 



(DIO) 
(DID 

(D12) 

(D13) 
(D14) 

(D15) 



If one takes the dot product of equation (3) with R/Bq, uses equation (4) to eliminate Bi and equation (D13) to eliminate pi, one 
finds that 

uirPo / _ mg\ ^ ma^poV-ui 

CO \ r J cor ' 



-iCOpoMie 



(D16) 



The term on the left-hand side represents inertia in the frame of reference rotating with the plasma, the terms Fqh = — Mir2i2po, 
Fg II = uirrngpo/cor, and FpU — — mfl^po(V • ui )/cor are the components of the Coriolis, gravitational, and pressure forces along the 
perturbed magnetic field line in roughly the +9 direction. One can show that for modes with complex CO, 






(D17) 



which is less than 1 for C, < ^^rit^ but which approaches 1 as ^ ^ ^^rit (as C, is increased towards ^c^t, a^G is increased towards g^). 
The stabilization of modes as ^ is increased to ^crit thus coincides with pressure becoming equal to gravity in the dynamics along 
field lines. With the use of equation (D12), equation (D16) can be rewritten as 



—iDsUiQ ~ —Mir I 2i2co ) 



(D18) 



Equations (D15) and (D18) show that the sign of Ds is critical for determining the phase of pi relative to uir- Let Fb be the sum of 
the components of the gravitational and Coriolis forces along the magnetic field. As explained in section 5, if CO is real and Ds < 0, 
the peaks in the perturbed density will occur at the points of convergence of Fb. If co is real and Dj > 0, the peaks in the perturbed 
density will occur at the points of divergence of Fb. 
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Fig. 1 . — A simplified model for an equilibrium of stt'ong vertical fields at the Galactic center. 
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Fig. 2. — The simplified one-dimensional equilibrium that is analyzed in this paper. 
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Table 1 
Notation 



Symbol Equation where introduced Meaning 

p, p, u, B (2), (3) mass density, flow velocity, pressure, and magnetic field; (0 

indicates equilibrium value, 1 indicates perturbed quantity) 

true gravity, effective gravity 

uniform angular velocity 

sound speed, Alfven speed, rotation speed 

adiabatic index 

azimuthal, axial, and radial mode numbers 

frequency in inertial frame, frequency in rotating frame 

convenient abbreviations 

dispersion-relation polynomials for the sound wave, 

Alfven wave, epicyclic-Alfven mode, and Rayleigh-Taylor mode 

ratio of thermal to magnetic pressure 

normalized radius 

exponent in constant-p equilibrium (po °^ r^'^) 

order-unity factor in constant-p equilibrium 

exponent of r in power-law solutions for u\r 

Lagrangian displacement 

measure of the wavenumber 

maximum value of ^ corresponding to an unstable mode 

perturbed gravitational, Coriolis, and pressure forces 

strength of the pervasive vertical field at the Galactic center 



g*,8 


(3), (11) 


Q. 


(8) 


a, VA, Vrot 


(13), (14), (78) 


Y 


(13) 


m, k^, n 


(15), (61) 


a, CO 


(15), (17) 


h, F, H, r, pn 


(23), (26), (42), (45), (63) 


D„ Da, D„ Dr 


(24), (27), (70), (71) 


P 


(30) 


X 


(34) 


q 


(33) 


a 


(39), (40) 


b 


(43) 


^ 


(47) 


c 


(73) 


Scrit 


(85) 


F„Fc,F, 


(86), (87), (D17) 


^vert 


(94) 
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Fig. 3. — Comparison of approximate analytic solution to dispersion relation for Rayleigh-Taylor modes [equation (80)] and numerical solutions of equations (62) 
and (68). The solid (dashed) lines represent the real (imaginary) part of co/fl from equation (80), while the crosses (diamonds) give the real (imaginary) part of 
numerical solutions of equations (62) and (68). The upper two panels correspond to equation (68), and the lower two panels correspond to equation (62). In the 
upper-right panel, q/va = 0.1 and va/vioi = 01> and the plasma is stable. In the other three panels, a/vA = 0.01 and va/vjo, = 0.2, and the plasma is unstable. For 
all of the plots, q = 'i. 
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Fig. 4. — Plots of DeZ)R (solid line) and D^Dr + S (dashed line), with S given by equation (74), a/i'A = 0.2, VA/vrot = 0.2, m = 10, n = 30, ln(r2/ri) = 5, and 
^ = 5. 
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Fig. 5. — The crosses indicate combinations of parameters (evaluated in the dense plasma) for which the numerical solutions of equations (62) and (68) indi- 
cate instability. The dots indicate parameters for which the numerical solutions are real. The solid line corresponds to the approximate analytic stability criterion 
equation (93). 
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density enhancement when 
Vphase > sound Speed, O) is real 



and IE II I > IE II I 
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Vphase < sound Speed, (O is real, 
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Fig. 6. — Eigenmodes corresponding to equation (80) with ni ~ o(l) and a Lagrangian total-pressure perturbation that vanishes everywhere. FgU and Fqu are 
the components of the gravitational and Coriolis forces along the perturbed magnetic field lines in the dense plasma. When (0 is real and Vp^ase < "J. the plasma is 
compressed at the points towards which Fg« +Fc.\\ converges, so that the perturbed pressure force opposes Fg« + Fqu. When CO is real and Vphasc > <^, the plasma 
responds inertially and is compressed at the points away from which Fg« + Fc,\\ diverges. In the latter case, magnetic tension is unable to balance the perturbed 
radial gravitational force, and no stable mode exists. 
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Fig. 7. — An m = 1 external kink mode. 
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Fig. 8. — Because the field lines in the hot low-density plasma are not anchored above or below the Galactic plane, they remain almost straight if they interact 
with a cloud at only a single point along their length provided the Alfven speed in the hot plasma v'A.hot is 2> Vdomj . 



